
#ifndef __LRWWSIMPLEX_STATE_H__

#include <gsl/gsl_vector.h>

/* defines a simplex in R^n */
typedef struct
{
  int n;
  gsl_vector **v;   /* (n+1)*n matrix of simplex vertices, 
                       one vertex per each matrix row */
  gsl_vector *fv;   /* (n+1)-dimensional vector of function 
                       values at the vertices */
  double V;         /* simplex volume */
}
simplex;

typedef struct
{
  simplex *S;
  gsl_vector *d;
  gsl_vector *xc;
  gsl_vector *xr;
  gsl_vector *xe;
  gsl_vector *xcc;
}
lrwwsimplex_state_t;

#define __LRWWSIMPLEX_STATE_H__

#endif
